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Abstract. The spin Seebeck effect refers to the generation of a spin voltage caused 
by a temperature gradient in a ferromagnet, which enables the thermal injection of spin 
currents from the ferromagnet into an attached nonmagnetic metal over a macroscopic 
scale of several millimeters. The inverse spin Hall effect converts the injected spin 
current into a transverse charge voltage, thereby producing electromotive force as in 
the conventional charge Seebeck device. Recent theoretical and experimental efforts 
have shown that the magnon and phonon degrees of freedom play crucial roles in the 
spin Seebeck effect. In this article, we present the theoretical basis for understanding 
the spin Seebeck effect and briefly discuss other thermal spin effects. 
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1. Introduction 

Generation of electromotive force by a temperature gradient has been known for many 
years as the Seebeck effect In recent years, a spin analogue of the Seebeck effect 
has drawn much attention in the field of spintronics, because replacing charge transport 
with spin transport in modern solid state devices is a major issue in the spintronics 
community. More than two decades ago, Johnson and Silsbee [2] published a seminal 
theoretical study, in which they generalized the interfacial thermoelectric effect to 
include spin transport phenomena. Because their framework implicitly relies on a spin 
transport carried by spin-polarized conduction electrons, the phenomenon discussed in 
Ref. [2] should be classified as a "spin-dependent" Seebeck effect from this perspective. 
The field of thermal spintronics is sometimes called spin caloritronics An experiment 
reported in 2008 put a new twist on spin caloritronics, because understanding of that 
experiment requires a framework other than the "spin-dependent" Seebeck effect. 

In 2008, Uchida et al. demonstrated that when a ferromagnetic film is placed under 
the influence of a temperature gradient, a spin current is injected from the ferromagnetic 
film into attached nonmagnetic metals with the signal observed over a macroscopic scale 
of several millimeters [I]. This phenomenon, termed the spin Seebeck effect, surprised 
the community because the length scale seen in the experiment was extraordinarily 
longer than the spin-flip diffusion length of conduction electrons, suggesting that the 
conduction electrons in the ferromagnet are irrelevant to the phenomenon. Subsequently, 
the spin Seebeck effect was observed in various materials ranging from the metallic 
ferromagnets Co2MnSi |5j to the semiconducting ferromagnet (Ga,Mn)As [6J, and even 
in the insulating magnets LaY 2 Fe 5 0i 2 [Z] and (Mn,Zn)Fe 2 4 [8]. These observations 
have established the spin Seebeck effect as a universal aspect of ferromagnets. 

In a spin Seebeck device, the spin current injected into an attached nonmagnetic 
metal is converted into a transverse charge voltage with the help of the inverse spin 
Hall effect j9j [101 [11]. Therefore, the spin Seebeck effect enables the generation of 
electromotive force from the temperature gradient as in conventional charge Seebeck 
devices. What is new in the spin Seebeck device is that it has a scalability different from 
that of conventional charge Seebeck devices, in that the output power is proportional 
to the length perpendicular to the temperature gradient. In addition, the paths of the 
heat current and charge current are separated in the spin Seebeck device in contrast 
to the charge Seeebck device, such that the spin Seebeck device could be a new way 
to enhance the thermoelectric efficiency. Because of these new features, an attempt is 
already underway to develop a new spin Seebeck thermoelectric device [12l HHJ [14] . 

As is inferred from the fact that the spin Seebeck effect occurs even in an insulating 
magnet [7] , this phenomenon cannot be described by the the "spin- dependent" Seebeck 
framework proposed by Johnson and Silsbee [2] . Instead, we need several new ideas and 
notions. In this article, we introduce basic ideas to understand the spin Seebeck effect. 
In addition, we present a brief summary of other thermo-spin phenomena. 
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2. Spin Current 

The spin Seebeck effect is a long-range thermal injection of the spin current from a 
ferromagnet into an attached nonmagnetic metal. Therefore, knowledge on the spin 
current is indispensable for understanding the spin Seebeck effect. In spin-orbit coupled 
systems, the spin is a non-conserved quantity, and hence there have been a number of 
discussions on the proper definition of spin currents in such systems [151 [16]. We do not 
discuss this subtle problem in this article, but here we present a simple argument. Let 
us consider the following definition of a spin current «7 S : 

Js = ^2s g k V k , (1) 

fc 

where s k is the z-component of the spin density s k with the z axis chosen as a spin- 
quantizing axis, and v k is the velocity of elementary excitations concomitant to the spin 
density s k . We consider here a spin-independent velocity v k because we focus on a pure 
spin current that is unaccompanied by a charge current. 

From Eq. (CQ) we can derive two kinds of pure spin current. The first is the so- 
called conduction-electron pure spin current. In this case, the ^-component of the 
spin density is given by s| = c k ^c k ^ — <4, iC*.^, where c ka is the creation operator for 
conduction electrons with spin projection o =f, \. and momentum k. After taking the 
statistical average, the expectation value of the conduction-electron pure spin current 
J s c_el is calculated to be 

J s" cl = E^((it c ^)-(4^))' ( 2 ) 
fc 

where v k is the velocity of conduction electrons. From this expression, we see that an 

asymmetry between the up-spin population and the down-spin population is necessary 

to obtain a nonzero conduction-electron pure spin current. 

The second type of pure spin current is the so-called magnon spin current. In this 

case the z-component of the spin density is given by s k — So — b k b k , where b k is the 

creation operator for magnons with momentum k. Substituting this into Eq. and 

taking the statistical average, the expectation value of the magnon pure spin current 
jmag j g gj ven 

^ = 4X>*(< 6 t 6 *> - (bl h b- k )), (3) 
fc 

where v k is the magnon velocity, and we have used the relation v_ k = —v k . From 
this expression, we see that an asymmetry between the left-moving population and the 
right-moving population is necessary to obtain a nonzero magnon spin current. 

These two spin currents can be detected experimentally in the following way. For 
the conduction-electron spin current «7g , the method of nonlocal spin injection and 
detection is used [TTJ, [18] . In the device shown in Fig. [TJ a charge current I c is applied 
across the interface between a metallic ferromagnet F\ and a nonmagnetic metal N. 
Because the conduction electrons in F\ are spin polarized, a spin-polarized current is 
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Figure 2. Schematic of a device that injects and detects the magnon spin current [20] . 



injected from F% into N, which creates a spin accumulation at the interface between 
Fx and N. Then, because there is no charge current flowing to the right side of Fi, 
the spin accumulation at the Fi/N interface diffuses to the right in the form of a 
conduction-electron spin current. The signal of the conduction-electron spin current is 
detected through the second metallic ferromagnet F 2 by measuring the electric voltage 
as shown in Fig. [TJ If and only if there is a spin accumulation at the F 2 /N interface, 
the electrochemical potential at the F 2 /N interface is influenced by whether or not the 
magnetization in F 2 is parallel to that in Fx (for more details, see Ref. [T5]). 

For a magnon spin current «J s mag , an insulating magnet is used to eliminate the 
contribution from the conduction-electron spin current [20J. In Fig. [2j two platinum 
films are put on top of a yttrium iron garnet (YIG) film. The first Pt film (Ni) acts as 
a spin current injector with the help of the spin Hall effect (see the next section). The 
spin current injected from Ni exerts spin torque on the localized magnetic moment at 
the Ni/F interface. Owing to the spin torque, the magnetization at the N\/F interface 
starts to precess and induces a spin current. Then the spin current propagates through 
F in the form of a magnon spin current «7 s mag . When the magnon spin current propagates 
from the N\/F interface to the N 2 /F interface and the localized spins at the interface are 
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Figure 3. (a) Definition of k^ ut = fc out — (fe in • fc out )fci n appearing in Eq. (O, where 
k[ n and /c ut are the incident and scattered wave vectors, (b) Schematic of the spin 
Hall effect. The charge current J c is converted into the transverse spin current J s . (c) 
Schematic illustration of the inverse spin Hall effect. The spin current J s is converted 
into the transverse charge current J c . The spin-quantizing axis is perpendicular to the 
plane of the sheet. 



excited, the spin current is injected from F into N2 owing to the s-d exchange interaction 
at the interface [21]. The spin current thus injected can be detected electrically via the 
inverse spin Hall effect (see the next section). 

The important point here is the difference in the decay lengths between the 
conduction-electron spin current and the magnon spin current. The conduction-electron 
spin current J s c_el decays over 10-100 nm in noble metals depending on the strength of the 
spin-orbit interaction [19]. On the other hand, the magnon spin current can sometimes 
propagate over a macroscopic length scale of a millimeter, which was indeed observed 
in Ref. [20]. 

3. Spin Hall Effect 

The spin Hall effect refers to the appearance of a nonzero spin current in the direction 
transverse to the applied charge current. In the spin Seebeck effect, the spin current 
injected from a ferromagnet into an attached nonmagnetic metal is converted into a 
transverse charge voltage via the reverse of the spin Hall effect, the so-called inverse 
spin Hall effect [91 [101 [TTj . Namely, the inverse spin Hall effect is used for electrical 
detection of the spin Seebeck effect. There are already a number of publications on the 
spin Hall effect in the literature, and we recommend Ref. [19] for readers interested in 
the detailed derivation of the spin Hall effect. 

The basic idea of the spin Hall effect [221 E3] is as follows. It is known that, in the 
presence of the spin-orbit interaction, a scattered electron acquires a spin polarization 
with the polarization vector a given by 

a oc k in x fc out , (4) 
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where fc in and fc out are the incident and scattered wave vectors. By multiplying both 
sides of Eq. by the vector k in , we see that the component of the scattered wave 
vector perpendicular to the incident wave vector, i.e., k^ ut = k out — (fc in • fc ou t)&in, is 
given by 

fc^ ut oc <x x k- m . (5) 

This equation means that the scattered vector is determined by the spin state and wave 
vector of the incident electrons. Macroscopically the spin Hall effect can be expressed 
as pEH El E51 EE] 

J s = 9 H a x J c , (6) 
while the inverse spin Hall effect is expressed as 

j c = e H a x j s , (7) 

where 9h is the spin Hall angle and a denotes the direction of the spin polarization. 

We now explain how the spin Hall effect works by taking the experiment of injection 
and detection of a magnon spin current via the spin Hall effect |20j (Fig. E]) as an 
example. Here, the nonmagnetic metal Ni is used as a spin-current injector by means 
of the spin Hall effect. In iVi, a charge current J c is applied parallel to the x direction. 
Then the spin current «7 S (|| y) across the Ni/F interface that is generated by the spin 
Hall effect has a spin polarization along the z axis owing to Eq. ([6]). This spin current 
«7 S creates a spin accumulation /x (|| z) at the Ni/F interface, and through the s-d 
exchange interaction at the interface [27J it exerts a spin torque on the magnetization 
M at the N x /F interface in the form TocMx(Mx/x) [2H EHl EI]- This torque 
excites a magnon spin current in the ferromagnet. 

The nonmagnetic metal N 2 is used to detect the magnon spin current by means of 
the inverse spin Hall effect. When the magnon spin current propagates from the N\/F 
interface to the N 2 /F interface, it injects spins from F into N 2 with a spin polarization 
along the z axis, again owing to the s-d exchange interaction at the interface [21]. The 
injected spins polarized parallel to the z axis diffuse along the y axis, and are converted 
into a charge current along the x axis owing to the inverse spin Hall effect Eq. ([7]). 
Therefore, the magnon spin current is detected as a charge voltage as shown in Fig. El 
The inverse spin Hall effect plays an important role in electrically detecting the spin 
Seebeck effect. 

4. Spin Seebeck Effect 

The spin Seebeck effect is the generation of a spin voltage caused by a temperature 
gradient in a ferromagnet. Here, the spin voltage is a potential for electrons' spin to 
drive spin currents. More concretely, when a nonmagnetic metal is attached on top 
of a material with a finite spin voltage, a nonzero spin injection is obtained. In this 
section, we first present a brief summary of the spin Seebeck effect, and then show the 
experimental details of this effect. 
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Figure 4. Schematic of the experimental setup for observing the spin Seebeck effect [3]. 
Inset: Spatial dependence of the observed voltage. The data are measured by changing 
the position of the Pt strip. 

4-1. Brief summary of the spin Seebeck effect 

Figure H] shows the experimental setup for observing the spin Seebeck effect in a magnetic 
insulator LaY 2 Fe50i2 [7j. Here a Pt strip is attached on top of a LaY 2 Fe50i2 film in 
a static magnetic field H = H z (3> anisotropy field), which aligns the localized 
magnetic moment along z. A temperature gradient VT is applied along the z-axis, 
which induces a spin voltage across the LaY 2 Fe50i2/Pt interface. Then this spin voltage 
injects a spin current I s into the Pt strip (or ejects it from the Pt strip). A part of the 
injected/ejected spin current I s is converted into a charge voltage through the inverse 
spin Hall effect 0Q11CEI]: 

V = 9 H (\e\I s )(p/w), (8) 

where \e\, 8h, p, and w are the absolute value of electron charge, the spin Hall angle, 
the electrical resistivity, and the width of the Pt strip (see Fig. @J. Hence, the observed 
charge voltage V is a measure of the injected/ejected spin current I s . By using this 
configuration, the spin Seebeck effect is observed not only in ferromagnetic metals (NiFe 
alloys [4] and Co2MnSi pj), but also in ferromagnetic semiconductors ((Ga,Mn)As) [6] 
and insulators (LaY 2 Fe 5 0i 2 [7] and (Mn,Zn)Fe 2 4 [8J). 

As shown in the inset of Fig. HI the spatial dependence of the spin Seebeck effect 
can be measured by changing the position of the Pt strip. Note that the signal has 
a quasi-linear spatial dependence, with the signal changing signs at both ends of the 
sample and vanishing at the center of the sample. 

It has been shown that the conduction electrons alone cannot explain the spin 
Seebeck effect, because the conduction electrons' short spin-flip diffusion length (~ 
several nanometers in a NiFe alloy) fails to explain the long length scale (~ several 
millimeters) observed in experiments (30JI3T]- This interpretation is further supported 
by the following two experiments. As we have already discussed, it was demonstrated 
in Ref. [20J using a ferromagnetic insulator YIG that spin currents can be carried 
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Figure 5. (a) AT dependence of V in the LaY2FesOi2/Pt sample at H = 100 Oe, 
measured when the Pt wires were attached near the lower-temperature (300 K) and 
higher-temperature (300 K+AT) ends of the LaY2Fe50i2 layer, (b) H dependence of 
V in the LaY 2 FesOi2/Pt sample for various values of AT. 



by magnon excitations. Subsequently, it was reported that, despite the absence of 
conduction electrons, the spin Seebeck effect can be observed in LaY 2 FesOi2, a magnetic 
insulator [7]. These experiments suggest that, contrary to the conventional view of the 
last two decades that the spin current is carried by conduction electrons [32], the magnon 
can be a carrier for the spin Seebeck effect. 

Now there is a consensus that the spin Seebeck effect is caused by a nonequilibrium 
between the magnon system in the ferromagnet and the conduction electron system in 
the nonmagnetic metal. In certain situations, both the nonequilibrium magnons and 
the nonequilibrium phonons play an important role. 

4-2. Experimental details of the spin Seebeck effect 

Here we show experimental data on the spin Seebeck effect in a LaY 2 FesOi2/Pt sample. 
The sample consists of a LaY2FesOi2 film with Pt wires attached to the top surface. 
A single- crystal LaY 2 Fe50i2 (HI) film with a thickness of 3.9 /xm was grown on a 
GdsGasO^ (HI) substrate by liquid phase epitaxy, where the surface of the LaY 2 FesOi2 
layer had an 8x4 mm 2 rectangular shape. Two (or more) 15-nm-thick Pt wires were 
then sputtered in an Ar atmosphere on the top of the LaY 2 Fe50i2 film. The length and 
width of the Pt wires were 4 mm and 0.1 mm, respectively. 

Figures 0(a) shows the voltage V between the ends of the Pt wires placed near 
the lower- and higher-temperature ends of the LaY2FesOi2 layer as a function of the 
temperature difference AT, measured when a magnetic field of H = 100 Oe was applied 
along the z direction. The magnitude of V is proportional to AT in both Pt wires. 
Notably, the sign of V for finite values of AT is clearly reversed between the lower- and 
higher-temperature ends of the sample. This sign reversal of V is characteristic behavior 



Theory of the Spin Seebeck Effect 



9 




z R (mm) 

Figure 6. Dependence of V on zp t , the displacement of the Pt wire from the center 
of the LaY2FesOi2 layer along the z direction, in the LaY2FesOi2/Pt sample at 
AT = 20 K. 

of the inverse spin Hall voltage induced by the spin Seebeck effect. 

As shown in Fig. EJb), the sign of V at each end of the sample is reversed by 
reversing H. It was also verified that the V signal vanishes when H is applied along 
the x direction, which is consistent with Eq. (j3J). This V signal disappears when the Pt 
wires are replaced by Cu wires with weak spin-orbit interaction. These results confirm 
that the V signal observed here is due to the spin Seebeck effect in the LaY 2 Fe 5 0i 2 /Pt 
samples. 

The spatial dependence of the thermally generated spin voltage can be measured 
by changing the position of the Pt strip. As shown in Fig. [61 the V signal induced by 
the spin Seebeck effect has a quasi-linear spatial dependence, with the signal changing 
signs at both ends of the sample and vanishing at the center of the sample. 

5. Linear-Response Theory of the Spin Seebeck effect 

5.1. Local picture of thermal spin injection by magnons 

As we have already discussed, the conduction electrons in the ferromagnet are considered 
to be irrelevant to the spin Seebeck effect. The fact that the spin Seebeck effect is 
observed even in a magnetic insulator suggests that the dynamics of localized spins in 
the ferromagnet, or magnon, is important to the spin Seebeck effect. To understand 
the spin Seebeck effect from this viewpoint, we first consider a model for the thermal 
spin injection by localized spins (see Fig. [7]). In this model we focus on a small region 
encircled by the dashed line, in which a ferromagnet (F) with a local temperature Tp and 
a nonmagnetic metal (N) with a local temperature T/v are interacting weakly through 
interface s-d exchange coupling J S( j. For simplicity we assume that the region in question 
(encircled by the dashed line) is sufficiently small such that the spatial variations of any 
physical quantities can be neglected, and that the size of the localized spin is unity. It 
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Figure 7. Side-view schematic of the system considered in Sec. 15.11 for thermal spin 
injection. 



is also assumed that each segment is initially in local thermal equilibrium; then, the s-d 
exchange interactions are switched on, and the nonequilibrium dynamics of the system 
is calculated. 

The physics of the ferromagnet F is described by the localized moment M, for 
which the dynamics is modeled by the Landau-Lifshitz-Gilbert equation: 



7 (if + h) - ^fs 



a 



x M H M x dtM, 



(9) 



where H = H z is the external field, 7 is the gyromagnetic ratio, a is the Gilbert 
damping constant, and M s is the saturation magnetization. In the above equation, the 
noise field h represents the thermal fluctuations in F. By the fluctuation-dissipation 
theorem [33, 34J, h is assumed to obey the following Gaussian ensemble [35] : 



and 



< W(t) > = 

< W{t)h v (t') > 



2ksTF<y 
ia\M s 



(10) 

(11) 



where a| = hj/M s is the cell volume of the ferromagnet. 

The physics of the nonmagnetic metal N is described by the itinerant spin density 
s, and its dynamics is modeled by the Bloch equation: 



d t s= 1 s 



M\ J sd M 

So — X S + i, 



(12) 



where r s f is the spin- flip relaxation time, and s = Xjv«^sd is th e local equilibrium spin 
density [27] with the paramagnetic susceptibility xn in N. In this equation, the noise 
source I is introduced [36] Gaussian ensemble 



and 



< p(t) > = 

< F(t)l v {t') > 



2k B T N XN 

T s f 



5fj,,u5(t - 1'), 



(13) 
(14) 



to satisfy the fluctuation-dissipation theorem [33l 134"] . From now on we focus on the 
spin-wave region, where the magnetization M fluctuates only weakly around the ground 
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state value M s z, and M/M s = z + m is established to separate small fluctuations m 
from the ground state value. 

The central quantity that characterizes the spin Seebeck effect is the spin current 
I s injected into the nonmagnetic metal N, since it is proportional to the experimentally 
detectable electric voltage via the inverse spin Hall effect [Eq. (jBJ)]. This quantity can 
be calculated as the rate of change of the itinerant spin density in N as I s =< d t s z {t) >. 
Performing the perturbative approach in Eq. ffT2|) in terms of J s d, we obtain 

hit) = < s + (t)m-(t') > t ^ t , (15) 

h 

where s^ 1 = s x ± is y and m ± = m x ± im y . Introducing the Fourier representation 
f(t) = J ^f u z~ lu>t and using the fact that the right hand side of Eq. ( !T5|) is only a 
function oft — t' in the steady state, we obtain 

h= J f |" ^ « S >=« »> (16) 

where the average <C ■ • • ^> is defined by < s+mj, >= 2ti5(uj + lo') <C s^mZ^ 

To evaluate the right hand side of Eq. ( fl6l) . the transverse components of Eqs. ([9]) 
and (fT2"j) are linearized with respect to s± and mj-. Then, to the lowest order in J s d, we 
obtain 

4 = . ! n f C + y'*. -rhi) (17) 



and 



m o, = ^ — ( 7^ + ■ Sd 



tu — to — laco \ —iu + r~ f 

where Uq = ^Hq, h = h x ±ih y , and I = l x ±il y . From the above equations, we see that 
s and m are affected by both the noise field h in F and the noise source I in N through 
the s-d exchange interaction J sc j at the interface. Substituting the above equations into 
Eq. ( fl6l) . the spin current injected into N can be expressed as 

I s = J" - J s back , (19) 
where /P um P and J back are respectively defined by 

jpump = _ J sdS f°° du u < jhfrfhZu > . 2Q . 
hr si J_ 00 2tt \co - u + iau\ 2 \iuj - r^ 1 ] 2 



.back _ _ ^Sd [°° duJ u ^ f " L " > 

J S ~~ fc? / 0_ I . . • ..111- . _-ll9- ^"U 



and 

7-back ^ 

ft 2 y_ 00 27r \u — ojq + iao;| 2 |ia; — r^ 1 ) 2 
We readily see in this expression that JP um P represents the spin current pumped into 
N by the thermal noise field h in F (the so-called pumping component [37] ) , while /back 
represents the spin current coming back into F from the thermal noise source / in N 
(the so-called backflow component [58]). Using the two fluctuation-dissipation relations 
[Eqs. (TTTj) and (|T3l) ]. the pumping and backflow components are finally calculated to be 
= _ ( GskB / h)TF (22) 
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and 

4 back = - (G s k B /h)T N , (23) 
such that the net contribution can be summarized in the single expression 

I S = -G S ^-(T F -T N ), (24) 
n 

where G s = _ 2aT rfW2i r°° ^ ( « \ w j2 , h and | M = ft • 

used. Here the negative sign before G s arises from defining the positive direction of I s . 
Equation (124")) is the basic equation to understand the spin Seebeck effect. 

At this stage it is important to note that, when the z component of the quantity 
< [m x dtm] z > is calculated from Eq. (J9j under the condition r^ 1 ^> coo and by 
neglecting the attachment of the nonmagnetic metal N, we can be show that the 
pumping component [Eq. (120]) ] can be expressed as 

jpump = < [ m x dtm y > _ ( 25 j 

On the other hand, from the above argument we observe that the backflow component 
is given by the same quantity evaluated at the local thermal equilibrium, i.e., 

1^ = -G s < [m x d t m] z > loc _ eq . (26) 

Therefore, the thermal spin injection by localized spins can alternatively be expressed 
as 

Is = ~G S ( < [m x d t m] z > - < [m x d t m] z >i oc _ cq ) . (27) 

This procedure was used in Ref. [39J to perform the numerical simulation of the spin 
Seebeck effect. 

Equations ffl9l) and ff24|) indicate that when both F and N are in local thermal 
equilibrium (i.e., Tp = Tjv) there is no net spin injection into the attached nonmagnetic 
metal N. However, conversely, this means that if the ferromagnet F deviates from 
the local thermal equilibrium for some reason, a finite spin current is injected into (or 
ejected from) the attached nonmagnetic metal N. This leads to the following simple 
picture for the spin Seebeck effect. Namely, the essence of the spin Seebeck effect is that 
the localized spins in the ferromagnet are excited by the heat current flowing through 
the ferromagnet, which then generates finite spin injections because of the imbalance 
between the pumping component JP um P and backflow component 7^ ack . It is important 
to note that the heat current that excites the localized spins has two contributions: the 
magnon heat current and the phonon heat current. Accordingly, there are two relevant 
processes in the spin Seebeck effect. The first, in which the localized spins are excited 
by the magnon heat current, corresponds to the magnon-driven spin Seebeck effect 
discussed in Refs. [ID] and [4T] . The second, in which the localized spins are excited by 
the phonon heat current, corresponds to the pho non-drag spin Seebeck effect discussed 
in Ref. HZI. 
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Figure 8. Schematic of system composed of ferromagnet (F) and nonmagnetic metals 
(N) divided into three temperature domains Fi/Ni, F 2 /N 2 , and F 3 ,7V 3 , with local 
temperatures T±, T 2 , T3, respectively. 



5.2. Linear-response approach to the magnon- driven spin Seebeck effect 

In the previous section we gave the local picture of thermal spin injection by magnons, 
but we did not discuss effects of magnon diffusion. Now, starting from the equivalent 
model to the previous subsection, we first reformulate the thermal spin injection by 
magnons in terms of the quantum many-body theory, and then extend it to thermal spin 
injection containing the effect of magnon diffusion jH]. Let us consider the model shown 
in Fig. |5J where the ferromagnet and the attached nonmagnetic metals are interacting 
weakly through interface s-d exchange coupling J s &. This model is essentially the same 
as that considered in the previous subsection (Fig. [7]). An important point in this model 
is that there is no local temperature difference between the ferromagnet and the attached 
nonmagnetic metals, i.e., = Tp l = T\, T/v 2 = Tp 2 = T 2 , and T/v 3 = Tp 3 = T 3 . It is 
assumed that each domain is initially in local thermal equilibrium without interactions 
with the neighboring domains. We then switched on the interactions between the 
domains, and calculate the nonequilibrium dynamics of the system. 

The localized spin in the ferromagnet is described by the exchange Hamiltonian 

(r it rj) ri 

where [r^Vj) denotes a pair of the nearest neighbors. In addition to Eq. (1251) . we 
later consider a selfenergy correction to represent the Gilbert damping term in the 
magnon propagator. The single-particle Hamiltonian for conduction electrons in the 
nonmagnetic metal is given by 

U N = ^2 c l{ e p S p,p' + Up-P'i 1 + i^socr ■ (p x p')]} c p'' (29) 
p,p> 

where c p = (Cp-pC^) is the electron creation operator for spin projection f and |, 
U p -pi is the Fourier transform of the impurity potential [/j mp Yl ro e impurities — r o); an d 
rj so is the strength of the spin-orbit interaction [19] . Finally at the F-N interface, the 
magnetic interaction between the conduction-electron spin density and localized spin is 
described by the s-d Hamiltonian, 

Kf-n = £ J* d q s k ■ S q , (30) 
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where s k = J2 P c ^ p +k (TC p ls the spin-density operator of conduction electrons, 
S q = ^7= J2 q S( r i) eiq r is the localized spin operator at the interface, and Np (Nn) 
is the number of lattice sites in F (N) in each domain. Here, J k A ~ q is the Fourier 
transform of J sd (r) = J sd Er eF-N interface a s S ( r ~ r o) with ^sd being the strength of the 
s-d exchange interaction. 

The spin current induced in the nonmagnetic metal iV, (i = 1,2,3) can be calculated 
as the rate of change of the spin accumulation in N iy i.e., I s (t) = ^2 reN . (d t s z (r, t)) = 
(dtS ko (t)) ko ^ , where (■ ■ •) means the statistical average at a given time t, and s k = 
y/N^Sk- Introducing the magnon operator 

s x (rd = yj^5>-« + a «y q ' r ^ ( 31 ) 



Q 



and 



(33) 

5 ° + ^E G S+Q eiQ ' ri ' ( 34 ) 



Q,Q 



the Heisenberg equation of motion for s z ko yields 

q,k 

where So is the size of the localized spins in F. Here have used the relation \s k , s^,} = 
±2s^ +k , and neglected a small correction term arising from the spin-orbit interaction, 
assuming that the spin-orbit interaction is weak enough in the vicinity of the interface. 
The statistical average of the above quantity gives the spin current 

I.® = y ~ 4 ^d q ^ > ft eC < a, t), (36) 

q.k 

where C kq (t,t') = —i(aq(t')s k (t)) measures the correlation between the magnon 
operator a+ and the spin-density operator s k = (s k — is y k )/2. Note that the time 
dependence of I s (t) vanishes in the steady state and hence is hereafter discarded. 
Introducing the frequency representation C k Jt—t') = f^C , ^ fc (w)e _lw ^*'\ adopting 

the representation [43J C = (o*|c*)> and usin S the relation C< = \[C K - C R + C A ), 
we obtain 

qk V2N F N N h 

for the thermal spin current I s in the steady state, where we have introduced the 
shorthand notation f = f°° 

Jul J — OO Z7T 

We first consider the process shown in Fig. [9] (a) where magnons travel around 
the ferromagnet F\ without sensing the temperature difference between F\ and F 2 . 



Y 7 sd Vb o ne / C?Ju) (37) 
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Figure 9. Feynman diagrams expressing the magnon-driven spin Seebeck effect [Hj. 
The thin solid lines with arrows (bold lines without arrows) represent electron 
propagators (magnon propagators). 



Using the standard rules of constructing the Feynman diagram in Keldysh space |44j, 
the corresponding interface Green's function Ck, q {w) for the correlation between the 
magnons in F% and the itinerant spin density in N\ can be written in the form 

= J %^ Xk(u;)X q (u), (38) 

where Nn and Np are the number of lattice sites in Nx and Fx- In the above equation, 
Xfe(w) is the spin-density propagator 

'Xfc(w),x*M 
o, x£{u 

while X q (u) is the bare magnon propagator 



*»M = ( „ " fi ' ) • (39) 



both of which satisfy the equilibrium condition: 

x£W = lx%(u)]*, X? (") = 2i bnxjf («) coth(^) (41) 

and 

= [X*{u)]\X«{u) = 2iImX»coth(^;). (42) 

The retarded component of Xfe( w ) is given by x& ( w ) = Xn/(1 + ^%k 2 — kur s f) [45] where 
Xtv, Aat, and r sf are the paramagnetic susceptibility, the spin diffusion length, and spin 
relaxation time. The retarded component of X q (u) is given by X q {oj) = (lu— Uq+iau)^ 1 , 
where a is the Gilbert damping constant and u q = r yH Q + D ex q 2 is the magnon frequency 
with D ex being the exchange stiffness. 
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Substituting Eq. (135]) into Eq. ( 157)) and using the equilibrium conditions [Eqs. (JUJ) 
and ( 142)) ]. we obtain the expression for the spin current injected into N\ 
T 4iV int J s 2 d So 



V2h 2 N N N F 
x coth( 



£ / Im X *(u;)ImXj(u;) 



q,k 
ftjjJ 



coth(- 



'2k B T Nl ' y 2k B T F /\ ' 

where iVjut is the number of localized spins at the interface. From Eq. (143)) . it is clear 
that no spin current is injected into the nonmagnetic metal Nx when Nx and Fx have 
the same temperature. 

The above result that the injected spin current vanishes when T Fl = Tn x originates 
from the equilibrium condition of the magnon propagator [Eq. ( 142)) ]. When the magnons 
deviate from local thermal equilibrium, the magnon propagator cannot be written in the 
equilibrium form, and it generates a new contribution. To see this, let us consider the 
process shown in Fig. [9] (b) where the magnons feel the temperature difference between 
F\ and F 2 through the following magnon interaction between Fx and F 2 : 

-^-Y, 2J ~* s *- s -*> (44) 



n 



F-F 



q,q' 



where J^. q ' is the Fourier transform of J ex {r) = J cx ^ 



ro£F-F interface 



a$5(r — r ). We 

now treat all of the magnon lines as a single magnon propagator 5X q (u) in the following 

way: 

~ ' ~ ~: q (uj)X q> (uj)X q (u). (45) 



F q' 

Then the propagator is decomposed into the local-equilibrium part and nonequilibrium 
part via jl6] 

(46) 



8X q {u) = 5X l q e \u) + 5X n q 



-eq. 



where 



5X l q ^(u) 



SX l - e i' R (u) 
0, 



SX l -^' K i 
5X l ~ e i'- 



(47) 



is the local-equilibrium propagator satisfying the local-equilibrium condition, i.e., 



5X l q eq ' A = [5X l q eq ' R ]* and 5X l q eq ' K = [5X l q eq > R - 5X l q ec >' A } coth(^) with 



while 



SX. 



n-eq. 



CO 



'0, 5X n - e ^ K (oo) 
0, 

is the nonequilibrium propagator with 5X q l ~ eq,K (uj) given by 



(4f 



(49) 



5X n-eq,K 



:*) = £ 



coth( 



2kftTp 2 ' 



coth( 



Too 



2k B T Fl 



)]■ 



(50) 
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When we substitute Eq. into Eq. (157)) and use Eq. (155)) with X q (co) replaced 
by SX q (oj), we obtain the following expression for the magnon-mediated thermal spin 
injection: 

T -4J s 2 d S , o(2J cx 5 ) 2 iV int iV^ 



x |X>)| 2 ImX*M[coth(^) - coth(^)], (51) 

where N( nt is the number of localized spins at the F1-F2 interface. The integration 
over u can be performed by picking up only magnon poles under the condition 
ahjj q <C k B T Fl , k B T Nl (which is always satisfied when magnon excitations are well 
defined), yielding fjmx* (o;)|X,f M| 2 ImX^(u;)[coth( 1 ^)-coth( 1 ^)] » ^5(«, q - 

u; g /)ImXfe (^q)[coth( 2A , B y )— coth( 2 , q T )]. With the classical approximation cotl^^^;) ps 
^S^, we obtain 

= AT int (jy ) XjV r sf (a/A jV ) 3 T _ 

where A is the size of F\ along the temperature gradient, and T2 = 

Jo dx !l d V [(i+a; 2 ) 2 +y 2 (2SoJ"exT S f/fe) 2 ] which is approximated as T 2 « 0.1426 (T 2 « 
0.337/i/2S'oJe X r s f) for 2S Q J cx r s{ /h~tl (for 2S J ex r s f/h ^> 1). Eq. (J52)) expresses the signal 
of the magnon-driven spin Seebeck effect. 

5.5. Length scale associated with the spin Seebeck effect 

We have already seen in Eq. (125)) that the pumping component is given by the quantity 
jpump _ ^Q s < [ m x <9 t m] z >. Using this result and the scenario of the magnon-driven 
spin Seebeck effect, let us calculate the spatial dependence of If ump (R) and discuss the 
length scale associated with the spin Seebeck effect. The starting point is the Landau- 
Lifshitz-Gilbert equation for a bulk ferromagnet written in the form [37] 

dtM(R,t) = - V • J M (R,t) 

+ M(R, t) x ( - 7 [Ho + *)] + Jf d t M ^ *)) > ( 53 ) 
where the M M component of the magnetization current J M is given by [48J 

J ^ = ^ [MxVjMf (54) 
with D ex being the exchange stiffness. Here the Greek indices refer to the components in 
spin space, and the Latin indices refer to the components in the real space. In Eq. (153)) 
the Gilbert damping factor a is an anisotropic tensor [49] to account for the difference 
between the transverse dynamics and longitudinal dynamics [50], and it is represented 
here as a = diag(a_i_, a±, a\\). Note that the transverse damping a± is relevant to the 
ferromagnetic resonance experiment, while information on the longitudinal damping a» 



< hURiWiRj,*) > = V 3 ,; M,1/ <WW - 0> (56) 
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is quite difficult to obtain from experiments. As before, the thermal noise field is given 
by the Gaussian white noise obeying 

< h»{Ri,t) > = (55) 

and 

2k B T(R i )a ^ ) 
ia%M 8 

where a^ u = a\\ for // = v = x,y and a^ u = a± for // = v = z [51]. We use 
again the spin- wave approximation M/M s = z + m and rotating- frame representation 
m^ 1 = m x ± im y . Using the transverse component of the Landau-Lifshitz-Gilbert 
equation f )53|) and taking its statistical average, the pumping current If nmp (R) = —G s < 
m x (R,t)d t m y (R,t) — m y (R,t)d t m x (R,t) > is calculated to be 

IT*(R)= -^(u;(-iV ri )+u;(-iV r2 )) 

x < m + (r 1 ,t)m~(r 2 ,t) > ri ,r 2 ^R (57) 

where uj(— iV) = ■yH Q + D cx (— iV) 2 . Because the above equation contains the gradient 
operator V r acting solely on one of the pairs in the correlator < m + (ri,t)m~ (r 2 ,t) >, 
it is useful to introduce the Wigner representation with respect to the spatial coordinate 
in the following manner 



R=-(n + r 2 ), r = r 1 -r 2 , (58) 

where R represents the center of mass coordinate, while r represents the relative 
coordinate. In this representation, we have the relation 

< m + { ri ,t)m-{r 2 ,t) > Pl , P2 _> H = -2 V < m z JR,t) > e i<?r , (59) 



where < m z q (R,t) >= —-^Y^k < m ^+K/2^ : ) m q-K/2^) > eiKR , an d we have 
introduced the Fourier transformation m~(r) = ^= ^ k m^e lk r . This allows us to 
represent the pumping current as 

ir mp (R) = 2G s J2"i < m z q (R,t) >, (60) 

where u q = ojq + D ex q 2 , and we have used the quasiclassical approximation \K\ <C \q\. 

To calculate the pumping current from Eq. fl60|) . we take the statistical average of 
the z component of the LLG equation ( |53l) : 

d t <m z (R,t) > = - V R - < J m *(R,t) > -2a|| J^Wq < m z q {R,t) > 

+ 3m < m + {R,t)-fh-(R,t) >, (61) 
where the last term is evaluated with the help of the Wigner representation (|58|) to give 

3m < m+{R,t)jh-{R,t) >= - 2a "^^ )7 , (62) 
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where we have used the Fourier representation in frequency space in the intermediate 
step of the calculation. 

We use the following assumptions to solve Eq. ( l61~i) in a closed form. First, we 
assume Fick's law of magnon diffusion, 

< J mZ (R, t) >= -VV R < m z (R, t) >, (63) 

where T> is the diffusion constant. Second, we introduce a wavenumber qo roughly 
corresponding to the thermal de Broglie wavenumber with kinetic energy k^T [54], 
which satisfies 

^2u q < m z q (R,t) >^ u qo < m z (R,t) > (64) 

where we have used < m z (R,t) >= J2 q < m Q (Rjt) >■ Substituting Eqs. (|63|) and (|64|) 
into Eq. ( IBTT) . we obtain 

(d t - DV^) < m z (R,t) >= -2a {l cu qo < m z (R,t) > - ^-f^ J7 , (65) 

where the right hand side represents the sink due to the longitudinal Gilbert damping 
(the first term) and source due to the heat bath (the second term). This equation can 
be solved in terms of the magnon distribution < m z (R,t) >. 

Now we evaluate the spatial dependence of the spin Seebeck effect. From Eq. ( !60l . 
thermal spin injection by localized spins is given by 

I s (R) = 2G s u qo ( < m z (R, t) > - < m z {R, t) > loc _ cq ) , (66) 

where we have considered the contribution from the backflow component [Eq. ( 1271) ] and 
used the approximation [Eq. f )64|) ]. Under the local equilibrium condition there is no 
magnon diffusion, and by setting the both sides of Eq. ( 165]) equal to zero, we calculate 
the local equilibrium magnon distribution to be 

< m z (R,t) > loc _ cq = -MM, (67) 

hw qo 

where we have used a 3 s M s = 7/1. This equation represents the classical limit of the 
magnon distribution function (e hlJJq o/ kBT — 1) _1 as it should because we neglect the 
quantum fluctuation in the fluctuation-dissipation relation (156|) . In a current-carrying 
steady state with magnon diffusion, we can set the time derivative equal to zero in 
Eq. (165|) . and by putting < m z (R,t) > — < m z (R,t) >ioc-cq=< Sm z (R,t) > we obtain 

< 5m z (R, t) >= 4" < Sm z (R, t) >, (68) 
-V 

where we have introduced a new length 

\ 2 m = V/(2a {l uj qo ). (69) 

As is clear from the fact that the thermal spin injection by localized spins is given 
by < 5m z (R,t) > [see Eq. (|66|) ]. A m corresponds to the length scale associated with the 
magnon-driven spin Seebeck effect. Physically, A m corresponds to the length associated 
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Figure 10. Schematic of the device structure used to detect the acoustic spin 
pumping [SSI I5H] ■ The dashed line represents the external phonon. The other lines 
mean the same as in Fig. |H1 

with magnon number conservation, or in other words it is an energy relaxation length 
for magnons. However, as we have already noted, experimental information on the 
longitudinal damping constant an is lacking, such that no reliable estimate of A m is 
available at the moment. In the case of the phonon-drag spin Seebeck effect, A m 
is replaced by A p , which corresponds to the length associated with phonon number 
conservation, or in other words it is an energy relaxation length of phonons. 

6. Phonon-Drag Contribution to the Spin Seebeck Effect 

Phonon drag is a well-established idea in thermoelectricity [551 156] . Back in 1946, 
in the context of thermoelectricity, Gurevich pointed out that thermopower can be 
generated by nonequilibrium phonons driven by a temperature gradient, which then 
drag electrons and cause their motions [57]. This idea, now known as phonon drag, 
has been established as the principal mechanism behind low-temperature enhancement 
of thermopower. Here, nonequilibrium phonons are the key. In this section, we first 
discuss acoustic spin pumping to understand the role of nonequilibrium phonons in 
the spin Seebeck effect. Then we present a microscopic approach to the phonon-drag 
contribution to the spin Seebeck effect. 

6.1. Acoustic spin pumping 

To understand the role of nonequilibrium phonons in the spin Seebeck effect, it is 
instructive to discuss so-called acoustic spin pumping [581 ES] because the spin Seebeck 
effect is a kind of thermal spin pumping. In the acoustic spin pumping experiment, 
a hybrid structure of a ferromagnet F and a nonmagnetic metal N are attached to 
a piezoelectric actuator that acts as a nonequilibrium-phonon generator (see Fig. [TO]) . 
When nonequilibrium phonons are generated from the piezoelectric actuator and interact 
with magnons in the ferromagnet, the magnons deviate from the equilibrium distribution 
through magnon-phonon interaction and inject spin current into the nonmagnetic metal. 
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We consider the interaction of exchange origin between magnons and phonons (the so- 
called volume magnetostrictive coupling [60J ) , since this has been shown to give the 
largest contribution [61]. The so-called single-ion magnetostriction [60] arising from the 
spin-orbit interaction [62] is assumed to be negligible, because if the latter coupling was 
relevant to the experiment in Refs. [58] and [59J, the resultant acoustic spin pumping 
should be seen at GHz frequencies instead of the MHz frequency at which the acoustic 
spin pumping is experimentally observed. However, we note that in the experiment of 
Ref. [63] . the single- ion magnetostriction [60] arising from the spin-orbit interaction [62J 
seems to be dominant. 

We start from the exchange Hamiltonian 

where J cx (Ri — -Rj) is the strength of the exchange coupling between the ions at Ri and 
Rj. The instantaneous position of the ion is written as Ri = rj + n(r, ; ) where the lattice 
displacement u(ri) is separated from the equilibrium position tv Up to the linear order 
in the displacement, the exchange Hamiltonian flTOl can be written in the form 

"Hex = ^ u; q a q a q + U 

mag— ph; 

(71) 

where co q = jHo + 2So J2s ^ex(<5) J2 q [1 — cos(q • 5)] is the magnon frequency with the 
lattice vector S = asS, and 

H mag _ ph = J2^/ a s) ■ Hn) - u( ri + 5)]5(r. J )5(r i + 6) (72) 

ri,S 

is the magnon-phonon interaction with the magnon-phonon coupling g given by 
VJ ex ((5) = (g/a s )S. 

In our case of acoustic spin pumping, the phonon is colored by a single wavenumber 
and frequency. The lattice displacement field u for a fixed wavenumber K is expressed 
as [61] u(ri,t) = iJ2K=±K ^KUK(t)e lK °' Vi , where the polarization vector is odd 
under the inversion K — > — K, and Uxif) can be expressed as UKif) = u K {t) +u^ K (t)* 
to satisfy UKif) = U-K(t)*. Note that the spatial average of [tt(rj)] 2 is given by 
([w(rj)] 2 ) av = 2\Uk \ 2 - Using this representation of the displacement vector and 
introducing the magnon operator a, a) [Eqs. ( 13T]) -( 1M|) ]. the magnon-phonon interaction 
becomes 

^2 ^K, q U K a\ +K a q , (73) 

q,K=±K 

where A K , 9 = gfw q (K ■ e K ) with g = J2s * ' Vi7 cx(<5)/E ( 5 ^ex(^)] being the 
dimensionless magnon-phonon coupling constant. Note that, up to the lowest order 
in u, only the longitudinal phonons couple to magnons when the phonons propagate 
along the symmetry axis of the crystal [61] . 

Now we consider the process shown in Fig. HOI in which nonequilibrium phonons 
interact with magnons and cause their nonequilibrium, thereby injecting a spin current 
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into the attached nonmagnetic metal. As before, when we treat the phonon-dressed 
magnon lines as a single magnon propagator 5X q (ou), it has the form 

5X q (u) = ^K, q \U K \ 2 X q {u)X q _ K {u-v K )X q {u), (74) 

K=±K 

where Vk = v p K is the phonon energy for the phonon velocity v p . When we substitute 
Eq. ( 174]) into Eq. f}37]) and use Eq. f )38|) with X q (u) replaced by d~X q (oj), we obtain the 
expression 

T V2h(J* d S ) ^ 



1 s 



N P N F /N int ^ 

k,q,K=±Ko 



A k , q (v K )AK, q \U K \ 2 (75) 



for the acoustic spin pumping, where the quantity Ak^ q {y) is defined by 

A*>)= / %™X^)ZmX^ K (u - v)\XK(uj)\ 2 
J uj 

^^-kf 25 - 00 ^^]" (76) 

which describes the correlation among the magnon, the phonon, and the itinerant spin 
density. Note that the acoustic spin pumping [Eq. (175]) ] is proportional to the square of 
the phonon amplitude |?7k| 2 - Therefore, the acoustic spin pumping is proportional to 
the power of the external sound wave. 

6.2. Phonon drag in the spin Seebeck effect 

In this subsection, we discuss the effect of nonequilibrium phonons on the spin Seebeck 
effect. In contrast to the previous subsection, the phonon in the this case is not an 
external field with a single color, but a statistical variable obeying Bose statistics. 
Therefore, it is necessary to represent the displacement field u with the phonon operator 
as 



«w = i£&Vra^ (6K + ^ yK '"- (77) 

where M ion is the ion mass and b^ K (Jok) is the phonon creation (annihilation) operator for 
wavevector K, is the polarization vector, and vk is the phonon frequency. Note that 
here and hereafter the polarization index ( is omitted, because we consider a situation 
where ( is not mixed with each other. Using this representation, the magnon-phonon 
interaction ( 172]) is expressed as 



7= y^?K, q B K a\ +K a q , (78) 
where Bk = bx + b J _ K is the phonon field operator, and the magnon-phonon vertex is 



"^mag-ph .... 

V q K 



given by T Kiq = 2S g Y,s \/ (S-e K )(5-K)[l-cos(q-6)] with the phonon velocity 



Now we discuss the phonon-drag contribution to the spin Seebeck effect. A natural 
guess is to replace \Uk\ 2 in Eq. ( 175]) with the deviation of the phonon distribution 
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Figure 11. Feynman diagrams representing the phonon-drag contribution to the spin 
Seebeck effect [42]. The dashed fine represents a phonon propagator. Meanings of the 
other lines are the same in Fig. [9] 



function from its local equilibrium value, namely, \Uk\ 2 — u p > — < n p >i oc _ eq . In 
the following we show that this captures the essence of the phonon-drag contribution to 
the spin Seebeck effect. For illustration, let us first consider the process shown in Fig.fTTI 
(a), where the magnons emit and absorb phonons while traveling around the domain F\, 
but neither the phonons nor magnons sense the temperature difference between Fx and 
i*2- The phonon-dressed magnon propagator 5X q (u) in Fig. [TT] (a) can be expressed as 

8X q (u) = X q (uj)£ q (u)X q (oj) (79) 

with the selfenergy due to phonons, 



2N ? A 

+ D A {v)f 1 X q _{uJ)+D K {v)X q _{u-)), (80) 

where f is the Pauli matrix in the Keldysh space, and we have introduced the shorthand 
notations w_ = to — u, q = q — K, and J = The bare phonon propagator in 

the above equation has the form 

where the retarded component is given by D^iy) = {y — ^x + i/Tp) -1 — (^ + ^K + i/Tp) _1 
with t p being the phonon lifetime [64] , and each component satisfies the equilibrium 
condition: 

D A K W) = [Dftiy)\\ D K K (v) = 2iImZ^)coth(^). (82) 

When the phonons are in thermal equilibrium, the phonon-dressed magnon 
propagator [Eq. ([T9~]) ] can be written in the local-equilibrium form 

5X q (u) = 6X l q e «(io) } (83) 
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where each component of the propagator satisfies the local-equilibrium condition, 
SX l q e i> A (u) = [5X l q e i> R (u)]* and 5X l q e ^ K (oo) = [5X l q e ^ R (uj) - 5X l q e i' A (oo)} coth(^), 
with the retarded component given by 

K F Jv 

x^( W .)+Df(^_(,.)}, (84) 

Using the same procedure to obtain Eq. (1451) . we calculate the injected spin current to 
be 



y/2h 2 N N N F J w ky ' qK ' 



q,k 

COttl ( oL rr ) - COth(- 



From this expression, we see that no spin current is injected into the nonmagnetic metal 
N\ when N\ and F\ have the same temperature. 

The above result that the injected spin current vanishes when Tp 1 = originates 
from the local-equilibrium condition of the magnons [Eq. ( 183]) ] which is derived from 
the equilibrium condition of the phonons [Eq. ( 182]) ]. When the phonons deviate from 
thermal equilibrium, the corresponding phonon propagator 5Dk{v) can be written in 
the form [46] 

6D K (u) = 5D l ^(u) + 5D n K e<1 (u), (86) 

where 5D l ^ q {i>) is the local-equilibrium propagator with local-equilibrium conditions 
SDfr^iv) = [5D l K e «> R (u)\* and 6D&> K {v) = [5D l ^ R (u) - 8D l ^ A {v)] coth(^), 
while 5D r £ eg (i>) describes the deviation from local equilibrium. When we allow such 
a nonequilibrium distribution of the phonons, the phonon-dressed magnon propagator 
cannot be expressed in the equilibrium form [Eq. (183]) ]. Instead, it is expressed in the 
form of Eq. ( 146]) with the nonequilibrium component 

\2 



K F Jv 



x 5D n K eq ' K (u) 



coth( — ) J — coth(- 



, . v ..7) 

This nonlocal propagator can give rise to a nontrivial contribution to the injected spin 
current. 

With the above in mind, let us next consider the phonon-drag process shown in 
Fig. [IT] (b), where the phonons sense the temperature difference between F\ and F2 
while the magnons do not. The phonon interaction between F\ and F2 is described 

by m 

p ^ +K 'Bk ■ B K ,, (88) 

F K K' 



H F-F 
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where Vtf +K ' is the Fourier transform of O p (r) = Q Y, ro eF-F interface - r o), and 

f2 = v/2~Kp/ 'Mi on with the elastic constant i^ p . The corresponding nonequilibrium 
phonon propagator 5Dk(v) is given by 

= |^ +K '| 2 5 K (z,) J D K ,(z,) J D K (^), (89) 



which can then be written in the form of Eq. (155)) : 5D l ^ q = ' 2l ^0i-f q %^ ') with 



1 2 



= E^I^Ml 2 ^^), (90) 



<*"•«, ~ c- ~r\"n-eq , K 

&nd5D n K eg = (° SD % )with 



AT' F 



x I^MHcothf— — ) - coth(— — )]. (91) 

When we substitute Eq. (J9~lj) into Eq. flHTj) and use Eqs. (1381) and fl37|) . we obtain 
the phonon-drag contribution to the injected spin current as 



F k,q,K,K' Jv 



x Im J D^( z/ )[ co th(-^) - coth(-^H], (92) 

where L = y/2(J^ d S )fllN int Nl n jN F , and A fe>q (i/) is defined in Eq. (j^g). After 
integrating over u; by picking up the magnon poles, Ah tq (y) is calculated to be 
A k M = 4>) + A®{u) with A®(u) = -I^IW* _ (i Xfc ( w ,)) 

and Agj,(z/) = -Ke^/ 2 ^ _ (jL Xfc ( Wfl )) ' (Sgt)-,^. Note that owing to 

the symmetry in the ^-integration, the leading term A^ q {y) does not contribute to 
the thermal spin injection. Then, we can perform the integration over v by picking 
up the phonon poles, j u | J D^(z/)| 2 Im J D^,(z/)[coth(^-) - coth(^)] = -itt p 5(v k - 

^')[coth(^-) - coth(^)], which yields 



2 



x A k , q (u K )[coth(^r-) - coth(-^)], (93) 

where Ujy = v p /as- 

The above expression, which is proportional to the phonon lifetime r p , gives the 
phonon-drag contribution to the spin Seebeck effect. After a rather lengthy calculation, 
Eq. (|93|) is transformed into 

^^-^(IjjJBrp, (94) 
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where the dimensionless constant r eff is given by T^ ff = ( j^r^ ) > the factor 
R = ° 2* x J *fa)°^A/as) measures ^ ne strength of the magnetic coupling at the F/N 
interface, and B = E\ ■ B 2 where E\ = ^-gpr Jo""^ du sinh "^ 2 ^ is a function of 
thermally-excited phonons with the Debye temperature Td = huo/k-Q, and £>2 = 
(t/Tm) 9 / 2 ^ t B TMT 5f ^3 jTm/t dv ^2_ is a f unc tion of thermally-excited magnons with T M 
being the characteristic temperature corresponding to the magnon high-energy cutoff. 

The important point of Eq. fl9^|) is that the spin Seebeck signal due to phonon 
drag is proportional to the phonon lifetime r p , because the carriers of the heat current 
in this process are phonons. Because the phonon lifetime is strongly enhanced at low 
temperatures (typically below 100 K) owing to a rapid suppression of the umklapp 
scattering, Eq. ( 1941) suggests that the spin Seebeck effect is enormously enhanced at 
low temperatures. In contrast, the signal at zero temperature should vanish because 
of the third law of thermodynamics. Therefore, the phonon-drag spin Seebeck effect 
must have a pronounced peak at low temperatures. Note that although the possibility 
of similar enhancement of the magnon lifetime in the magnon-driven spin Seebeck effect 
[Eq. fl52l) ] is not definitely excluded, judging from the ferromagnetic resonance linewidth 
in Y 3 Fe50i2 [66J as a measure of the inverse magnon lifetime, it does not seem likely. 

To date, there are two experimental findings that support the existence of the 
phonon-drag spin Seebeck effect. The first is the observation of the predicted low- 
temperature peak in the temperature dependence of the spin Seebeck effect [63, [69] . In 
Ref. [12] the earliest experimental data on the spin Seebeck effect in LaY 2 Fe50i2 were 
theoretically analyzed, and the theory predicted that the spin Seebeck effect must show a 
pronounced peak at low temperatures as is discussed above. In Ref. [67] the temperature 
dependence of the spin Seebeck effect was measured in (Ga,Mn)As, and the data showed 
a pronounced peak at low temperatures consistent with the theoretical prediction [42J. 
In Ref. [69J the same trend was confirmed for YIG. The other experimental finding that 
supports the phonon-drag spin Seebeck effect is the observation of a spin Seebeck effect 
that is unaccompanied by a global spin current. Reference [6] reported that cutting 
the magnetic coupling in (Ga,Mn)As while maintaining the thermal contact allowed the 
spin Seebeck effect to be observed even in the absence of global spin current flowing 
through (Ga,Mn)As. The phonon-drag spin Seebeck effect can explain the "scratch" 
test experiment [6], although the idea of a magnon-driven spin Seebeck fails to explain 
the experiment. Moreover, in a recent study [58], an isolated NiFe alloy on top of a 
sapphire substrate was used to measure the spin Seebeck effect. This study excluded 
the possibility of a dipole-magnon-driven spin Seebeck effect for the "scratch" test 
experiment [6], and found that only the phonon drag by the substrate phonons could 
explain the experiment. One important point is that the experiment of Ref. [HE] was 
performed at room temperature; nevertheless, the spin Seebeck effect was observed 
with the signal extended over several millimeters, as in the first observation of the spin 
Seebeck effect in NiFe alloy [3]. This shows that the phonon drag can contribute to the 
spin Seebeck effect even at room temperature. 
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Figure 12. Schematic of the experimental setup for (a) the transverse spin Seebeck 
effect and (b) the longitudinal spin Seebeck effect. 



7. Varieties of the spin Seebeck effect 

7.1. Longitudinal spin Seebeck effect 

Up to this point, we have discussed the transverse spin Seebeck effect [Fig. [12] (a)], 
in which the direction of the thermal spin injection into a nonmagnetic metal is 
perpendicular to the temperature gradient. There is another type of spin Seebeck 
effect called the longitudinal spin Seebeck effect [HI EE] EH] [Fig. [12] (b)], in which 
the direction of the thermal spin injection into a nonmagnetic metal is parallel to 
the temperature gradient. While both conducting and insulating ferromagnets can 
be used for the transverse spin Seebeck effect, the longitudinal spin Seebeck effect is 
well defined only for an insulating ferromagnet because of the parasitic contribution 
from the anomalous Nernst effect [5] [70] [71]. The longitudinal spin Seebeck effect 
has been observed in mono crystalline [68] and polycrystalline [69] YIG (Y 3 Fe50i2) 
as well as in polycrystalline ferrite (Mn,Zn)Fe204 [8]. The longitudinal spin Seebeck 
effect is the simplest configuration in which a bulk polycrystalline ferromagnet can be 
used. Therefore, it is considered to be a prototype of the spin Seebeck effect from an 
application viewpoint. 

In Fig. [131 we show typical experimental results for the longitudinal spin Seebeck 
effect. The sample consists of a mono crystalline YIG slab and a Pt film attached to a 
well-polished YIG (100) surface. The length, width, and thickness of the YIG slab are 6 
mm, 2 mm, and 1 mm, while the corresponding dimensions of the Pt film are 6 mm, 0.5 
mm, and 15 nm. An external magnetic field H (with the magnitude H) was applied in 
the x-y plane at an angle 6 to the y direction (see Fig. [TBTa)). A temperature difference 
AT was applied between the top and bottom surfaces of the YIG/Pt sample. Figure 
[TBTa) shows the voltage V between the ends of the Pt layer in the YIG/Pt sample as 
a function of AT at H — 1 kOe. When H was applied along the x direction [9 = 90°), 
the magnitude of V was observed to be proportional to AT. The sign of the V signal 
at finite values of AT is clearly reversed by reversing the VT direction. The V signal 
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AT(K) H(kOe) 

Figure 13. (a) AT dependence of V in the YIG/Pt sample at H = 1 kOe, measured 
when VT was applied along the +z and — z direction. The magnetic field H was 
applied along the x direction (9 = 90°) and the y direction (8 = 0). (b) H dependence 
of V in the YIG/Pt sample for various values of AT at 9 — 90°, measured when VT 
was along the +z direction. 

also changes its sign with reversing H when 9 = 90° [Fig. [131(b)] and disappears when 
H is along the y direction (6 = 0) [Fig. [T3T a)]. These results are consistent with the 
symmetry of the inverse spin Hall effect induced by the longitudinal spin Seebeck effect 
[see Eqs. © and ©]. 

A major feature of the longitudinal spin Seebeck effect is that the sign of the spin 
injection is opposite to that in the transverse spin Seebeck effect, as shown in Fig. [T2j 
Focusing on the spin current injected into the nonmagnetic metal (N) close to the cold 
reservoir, the pumping component /P um P is less than the backflow component /^ ack in the 
transverse spin Seebeck effect. In contrast, JP um P is greater than J^ ack j n the longitudinal 
spin Seebeck effect. 

A linear-response approach to the longitudinal spin Seebeck effect was developed 
in Ref. (72]. Here, we present a phenomeno logical argument. First, recall that the spin 
Seebeck effect can be understood in terms of the imbalance between the thermal noise 
of the magnons in the ferromagnet and the thermal noise of the conduction-electron 
spin density in the nonmagnetic metal. The former noise injects the spin current into 
the nonmagnetic metal, while the latter ejects the spin current from the nonmagnetic 
metal. Because the thermal noise in each element can be related to its effective 
temperature through the fluctuation-dissipation theorem, the spin Seebeck effect can 
also be interpreted in terms of the imbalance between the effective temperature of the 
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magnons in the ferromagnet and the effective temperature of the conduction-electron 
spins in the nonmagnetic metal [see Eq. 024p] . 

Then, the signal sign reversal between the longitudinal and the conventional 
transverse spin Seebeck effects may be explained by the following conditions: (i) Most 
of the heat current in the ferromagnet/nonmagnetic-metal hybrid system at room 
temperature is carried by phonons (see discussion in Ref. [73] in the case of YIG), 
and (ii) the interaction between the phonons and the conduction-electron spins in the 
nonmagnetic metal N is much stronger than the magnon-phonon interaction in the 
ferromagnet F. In the longitudinal spin Seebeck experiment, the nonmagnetic metal is 
in direct contact with the heat bath, and thereby is exposed to the phonon heat current 
due to condition (i). Then, because of condition (ii), the conduction-electron spins in 
the nonmagnetic metal N are heated up faster than the magnons in the ferromagnet 
F, and the effective temperature of the conduction-electron spins in the nonmagnetic 
metal rises above that of the magnons in the ferromagnet F. In the conventional spin 
Seebeck setup, by contrast, the nonmagnetic metal N is out of contact with the heat 
bath, and the phonon heat current does not flow through the nonmagnetic metal N, 
while the ferromagnet F is in contact with the heat bath. This results in an increase 
in the effective magnon temperature in the ferromagnet F. Therefore, in this case, the 
effective temperature of the conduction-electron spins in the nonmagnetic metal N is 
lower than that of the magnons in the ferromagnet F. This difference can explain the 
sign reversal of the spin Seebeck effect signal between the longitudinal and transverse 
setups. 

7.2. Thermoelectric coating based on the spin Seebeck effect 

The spin Seebeck effect in magnetic insulators can be used directly to design thermo-spin 
generators and, combined with the inverse spin Hall effect, thermoelectric generators, 
allowing new ways to improve thermoelectric generation efficiency In general, the 
efficiency is improved by suppressing the energy lossess due to heat conduction and 
Joule dissipation, which are realized respectively by reducing the thermal conductivity 
k for the sample part where heat currents flow and by reducing the electrical resistivity 
p for the part where charge currents flow. However, the Wiedemann-Franz law 
(n e p = constant) limits this improvement in electric conductors when k is dominated 
by the electronic thermal conductivity K e . In contrast, k and p are independent in 
the spin Seebeck device since the heat and charge currents flow in different parts of 
the sample: k is the thermal conductivity of the magnetic insulator, and p is the 
electrical resistivity of the metallic wire. The spin Seebeck effect in insulators allows us 
to construct thermoelectric devices operated by an entirely new principle, although the 
thermoelectric conversion efficiency is still small at present. 

In 2012, Kirihara et al. proposed a new thermoelectric technology based on the spin 
Seebeck effect called "spin-thermoelectric (STE) coating" [12], which is characterized 
by a simple film structure, convenient scaling capability, and easy fabrication (Fig. 
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Figure 14. Concept of the STE coating based on the spin Seebeck effect [Hj. The 
STE coating exhibits a straightforward scaling: a larger film area leads to a larger 
thermoelectric output. Such a simple him structure can be directly coated onto heat 
sources with different shaped (curved or uneven) surfaces. 



IT4l) . In their experiments, an STE coating with a 60-nm-thick Bi-substituted YIG film 
was applied by using metal organic decomposition on a nonmagnetic substrate. Notably, 
thermoelectric conversion driven by the longitudinal spin Seebeck effect was successfully 
demonstrated under a temperature gradient perpendicular to such an ultrathin STE- 
coating layer (amounting to only 0.01 % of the total sample thickness). The STE 
coating was found to be applicable even to glass surfaces with amorphous structures. 
Such a versatile implementation of thermoelectric function may give rise to other ways 
of making full use of omnipresent heat. 

7.3. Position sensing via the spin Seebeck effect 

The longitudinal spin Seebeck effect in magnetic insulators has also been used in two- 
dimensional position sensing using a YIG slab covered with a Pt-film mesh [13] . Figure 
[T5l shows a schematic of the YIG-slab/Pt-mesh sample. When part of the sample 
surface was heated, the position of the heated part was found from the spatial profile of 
the spin Seebeck voltage in the Pt mesh. The advantages of two-dimensional position 
sensing using the spin Seebeck effect are the simplicity of the device structure and the 
production cost; this device structure can be made simply by fabricating a patterned film 
on a commonly-used sintered polycrystalline insulator. Therefore, this position-sensing 
method gives us a realistic application of the spin Seebeck effect in thermally-driven 
user-interface devices and image-information sensors. 
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Figure 15. A schematic of the YIG-slab/Pt-mesh structure. In Ref. [T3], part of 
the sample was heated by laser light, and the two-dimensional position information 
of the heated part was found by calculating the tensor product of the spatial profiles 
of the SSE voltage along the x and y directions. Here, an external magnetic field 
was applied along the diagonal (45°) direction of the Pt mesh for generating the spin 
Seebeck voltage in both directions. 



8. Other thermal spintronic effects 

So far we have focused on the spin Seebeck effect. Besides the spin Seebeck effect, there 
are several intriguing phenomena in which the interplay of spin and heat plays a crucial 
role. In this section we briefly review other thermal spintronic effects. 

8.1. Spin injection due to the spin- dependent Seebeck effect 

A thermally driven pure spin-current injection across a charge-conducting interface 
has recently been reported by several groups, in which the "spin-dependent Seebeck 
effect" plays an important role. Slachter et al. [H] demonstrated thermally driven pure 
spin-current injection and its electrical detection using the nonlocal lateral geometry 
of NiFe/Cu. The physics behind this experiment is based on the spin-dependent 
thermoelectric effect. The spin-dependent current j^^ is described by 



where a^, fi^, and are the spin- dependent conductivity, spin-dependent 

electrochemical potential, and spin- dependent Seebeck coefficient, respectively. The 
spatial distribution of the spin accumulation m — fi± is described by the Valet-Fert spin 
diffusion equation: 



where A is the spin-flip diffusion length. The essence of the experiment can be seen by 
solving these two equations under an appropriate temperature distribution across the 
NiFe/Cu interface. 

Le Breton et al. [75] demonstrated thermal spin injection from NiFe into Si through 
an insulating tunnel barrier SiC^/A^Os and called the phenomenon "Seebeck spin 
tunnelling". Here the injected spin current was detected by the Hanle effect, and the 




(95) 



(96) 



Theory of the Spin Seebeck Effect 



32 



observed signal was analyzed in terms of the "spin-dependent Seebeck effect". It is 
important to note that the direction of the spin injection in these two experiments is 
parallel to the temperature gradient, such that the signal could contain the contribution 
from the longitudinal spin Seebeck effect. 

From the Kelvin relation H^ = TS^^ with the spin-dependent Peltier coefficient 
II^j,, we expect the reciprocal process, i.e., the spin- dependent Peltier effect. Flipse et 
al. [76] have recently reported observation of this effect. 

8.2. Seebeck effect in magnetic tunnel junctions 

Several groups have measured the tunneling magneto-thermopower ratio of magnetic 
tunnel junctions, which was discussed analytically [77] and computed by a first- 
principles calculation [78J. Walter et al. [79] observed the tunneling magneto- 
thermopower in a CoFe/MgO/CoFe magnetic tunnel junction, and Liebing et al. [80] 
in a CoFeB/MgO/CoFeB junction. The signal is caused by the spin-dependent Seebeck 
effect. 

8.3. Magnon-drag thermopile 

It is well known that two drag effects contribute to the thermoelectric effect in magnetic 
metals: one is the phonon drag in which nonequilibrium phonons transfer momentum 
to conduction electrons to produce thermopower, and the other is the magnon drag 
in which nonequilibrium magnons transfer momentum to conduction electrons |55j. 
However, the magnon-drag effect is easily masked by the phonon-drag effect, and in 
general, it is quite difficult to investigate only the magnon-drag effect. Costache et 
al. [81] recently overcame this difficulty and proposed a device named the "magnon- 
drag thermopile" which provides information about the magnon-drag effect. The device 
is composed of many pairs of NiFe wires connected electrically in series with Ag wires, 
but placed thermally in parallel. When the two magnetizations in a pair of NiFe wires 
are in the parallel configuration, the thermopower is zero because the contributions 
of each wire are of the same magnitude but opposite signs. However, when the two 
magnetizations in a pair of NiFe wires are in the antiparallel configuration, there is a 
difference in the magnon states between the two wires, and the resultant thermopower 
is nonzero. Note that, in principle, although any electron-magnon scattering process 
other than the magnon drag can contribute to the observed thermopower, the magnon 
drag can dominate the signal when the energy dependence of the electron lifetime is 
negligible. 

8.4- Thermal spin-transfer torque 

Thermal spin-transfer torque is also a highly debated topic. Hatami et al. [82] discussed 
the thermal spin-transfer torque in magnetic nanostructures of metals, and Jia et al. [83] 
recently developed a first-principles estimation of the same process. This effect is 



Theory of the Spin Seebeck Effect 



33 



relevant to the thermally driven domain wall motion discussed analytically by Kovalev et 
al. [51] and computed numerically by Yuan et al. [85J Thermal spin-transfer torque has 
also been discussed in the context of magnetic insulators. Slonczewski [86] discussed 
the thermal spin-transfer torque resulting from the longitudinal spin Seebeck effect 
in ferrite. Spin-transfer torque caused by magnons is called magnonic spin-transfer 
torque [87] , and Hinzke et al. |88j discussed the role of thermal magnonic spin-transfer 
torque. Experimentally, an evidence for the thermal spin-transfer torque was reported 
by Yu et al. [89]. 

8. 5. Effects of heat current on magnon dynamics 

Another interesting subject is the dynamics of magnon wavepackets under the influence 
of a temperature gradient. Padron- Hernandez et al. [90] found that magnon wavepackets 
propagating along a YIG film are amplified when a temperature gradient is applied 
perpendicular to the YIG film. This experiment implies that the magnon damping term 
is canceled by the action of the temperature gradient, which leads to an amplification of 
the magnon wavepacket. The observed result was interpreted by the authors in terms 
of the magnonic spin-transfer torque of thermal origin in the longitudinal spin Seebeck 
configuration. 

Lu et al. [91] studied the effects of heat current on ferromagnetic resonance. Using a 
trilayered structure consisting of a micron-thick YIG film grown on a submillimeter-thick 
gadolinium gallium garnet substrate and capped with a nanometer-thick platinum layer, 
they found that a temperature gradient over the trilayer can control the ferromagnetic 
relaxation in the YIG film. The result was interpreted by the authors in terms of the 
magnonic spin-transfer torque of thermal origin. 

8.6. Anomalous Nernst effect and spin Nernst effect 

The anomalous Nernst effect refers to the generation of a voltage gradient || mxVT 
by applying a temperature gradient VT in a ferromagnetic material with a magnetic 
polarization vector m. This phenomenon has been studied systematically in various 
ferromagnetic metals by Miyasato et al. [92J, in (Ga,Mn)As by Pu et al. |93j, and in 
NiFe lateral spin valve by Slachter et al. |94j. It is important to note that if there is a 
thermal conductivity mismatch between the substrate and the ferromagnetic film when 
measuring the transverse spin Seebeck effect for a conducting magnet, there can be a 
parasitic contribution from the anomalous Nernst effect as pointed out in Ref . [5] . This 
issue was recently discussed again in Ref. [70J. 

The spin Nernst effect refers to the generation of a transverse spin current J s with 
the spin polarization a by a temperature gradient, i.e., J s || er x VT. Reference [95] 
theoretically discusses the the spin Nernst effect in a two-dimensional Rashba spin-orbit 
system under a magnetic field, and Refs. [96J, [97] the same effect in a zero magnetic field. 
The spin Nernst effect of extrinsic origin is analyzed through first-principle calculations 
in Ref. [98]. 
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8. 7. Thermal Hall effect of phonons and magnons 

When the time-reversal symmetry is broken by a magnetic field or magnetic ordering, a 
finite Hall response can occur in principle even in the case of charge-neutral excitations 
such as phonons and magnons. Recently, the thermal Hall effect of phonons and 
magnons has been reported. Strohm et al. observed the thermal Hall effect of phonons 
in a paramagnetic insulator of terbium gallium garnet [99]. The result was explained 
by the interaction of local magnetic ions with the local orbital angular momentum of 
oscillating surrounding ions |100[ 1101] . The thermal Hall effect of magnons is also 
observed in an insulating ferromagnet LU2V2O7 with pyrochlore structure |102j . and the 
result was explained in terms of a Dzyaloshinskii-Moriya interaction. The Hall effect of 
magnons was also discussed theoretically in Refs. |103[ I104| [105] . 

9. Conclusions and future prospects 

We have discussed the physics of the spin Seebeck effect and clarified the important 
role played by magnons. Moreover, we have shown that nonequilibrium phonons also 
play an active role. From a theoretical viewpoint, one of the open questions in the 
spin Seebeck effect is the role of spin-polarized conduction electrons in the metallic 
and semiconducting ferromagnets, especially in interpreting of the experiment reported 
in Ref. |106] . From an experimental viewpoint, on the other hand, one of the most 
challenging issues is clarifying to what extent the spin Seebeck effect can be applied. A 
small but a firm step is already in progress [121 EES EH] • 

More generally, one of the driving forces for investigating thermal effects in 
spintronics is the desire to deal with heating problems in modern solid-state devices. 
From this viewpoint, thermo-spintronics is still in its infancy, and many issues still 
remain unclear. For example, the relationship between the pure spin current and 
dissipation |107] needs to be investigated extensively. Although the practical application 
of thermo-spintronics looks remote at present, we can definitely say that the interplay 
of spin and heat manifests itself in state-of-the-art experiments and involves interesting 
physics. 
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